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Distinguishing assigned quantum states with assigned probabilities via quantum measurements 
is a crucial problem for the transmission of classical information through quantum channels. Mea¬ 
surement operators maximizing the probability of correct discrimination have been characterized by 
Helstrom, Holevo and Yuen since 1970’s. On the other hand, closed-form solutions are available 
only for particular situations enjoying high degrees of symmetry. As a suboptimal solution to the 
problem, measurement operators, directly determined from states and probabilities and known as 
square root measurements (SRM), were introduced by Hausladen and Wootters. These operators 
were also recognized to be optimal for pure states equipped with geometrical uniform symmetry 
(GUS). In this paper we discuss the optimality of the SRM and find necessary and sufficient con¬ 
ditions in order that SRM maximize the correct decision probabilities for set of states formed by 
several constellations of GUS states. The results are applied to some specific examples concerning 
double constellations of quantum phase shift keying (PSK) and pulse position modulation (PPM) 
states, with possible applications to practical systems of quantum communications. 

PACS numbers: 03.67.Hk 


I. INTRODUCTION 

The long standing problem of transmitting classical 
information through a noiseless quantum channel refers 
to the following scenario. A classical information source 
emits a symbol S in a finite alphabet, say {sq, • ■ ■, Sm-i}, 
whose nature is irrelevant to the problem. The symbols s, 
are emitted with probabilities qi satisfying the properties 
gi > 0 and ft = 1- the basis of the symbol Si 

emitted by the classical source, the transmitter (conven¬ 
tionally known as Alice) prepares the quantum channel 
in a state belonging to a set of agreed quantum states, 
generally mixed states with density operators pi which 
are semidefinite positive {pi > 0) and have unitary trace 
(TrP j = 1). The quantum ensemble {qi,Pi\ can be con¬ 
veniently represented by the weighted density operators 
Wi = qiPi- 

The problem of the receiver (conventionally known as 
Bob) is to extract as well as he can classical information 
from the quantum channel he and Alice share. In this 
task Bob is assumed to know both the quantum states 
Pi and their probabilities g^. We assume that the mea¬ 
surement and decision strategy of Bob is the quantum 
hypothesis testing introduced by Helstrom [T]. The strat¬ 
egy consists in performing on the quantum channel a pos¬ 
itive operator-valued measurement (POVM) with mea¬ 
surement operators Pi semidefinite positive {Pi > 0) and 
resolving the identity operator, namely, ~ 

The conditional probability p {j\i) that the measurement 
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gives the result R = j provided that the state is pi turns 
out to be Tr (piPj). The usual task of Bob is the choice 
of the POVM Pi maximizing the probability of correct 
decision 

m — 1 m—1 m—1 

Pc = X ftp(*l*) = X ftTr(piPz) = X Tr(lTiPj) . 
2=0 2=0 2=0 

( 1 ) 

The problem solution is given by well known general re¬ 
sults due to Holevo [5], Yuen et al. [3], and Helstrom [T] 
(see also jl], and for a connection with min-entropy 
of a classical source) and is summarized by the following 
theorem. 

Theorem 1. The measurement operators Pi give the 
maximum correct detection probability if and only if the 
operator 

m — 1 

y = ^ w,p, (2) 

i =0 

satisfies the conditions 

y - Wi > 0 (3) 

for each f = 0,..., m — 1. The above condition implies 
also that Y is semidefinite positive and that the equalities 

P,{W.-W,)P,=0 ( 4 ) 

hold true for alH, j = 0,..., m — 1. □ 

Unfortunately, closed form solutions are available only 
in few very particular cases, generally implying pure 
states enjoying high degrees of symmetry. Indeed, as 
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pointed out in |^, serious mathematical difficulties arise 
because of the non linearity of the problem, which turns 
out to be a convex semidefinite programming problem. 

If the states prepared by Alice are pure, namely, pi = 
l7i)(7i|, it can be shown [7] that also the measurement 
operators have rank 1, namely. Pi = \p,i){p,i\. Then the 
correct detection probability Q depends on the inner 
products of the states and of the measurement vectors 
\pi), namely, 

m— 1 

(5) 

i=0 

where ly^) = y^|yi) are the weighted states. The opti¬ 
mization problem reduces to find the measurement vec¬ 
tors \p,i) maximizing the correct detection probability 
([^. In accordance with the general case, the vectors 
\pii) must resolve the identity operator I of the Hilbert 
space spanned by the states 

m— 1 

'^\pi){pi\= I. (6) 

i=0 

In the following we assume that the weighted states |yi) 
are linearly independent. This is not at all too restrictive 
for all situations of practical interest. In this case it can 
be shown |S] that the optimal measurement vectors \pi) 
(to be determined) form an orthonormal basis. Since the 
weighted states are linear combinations of the measure¬ 
ment vectors, provided that the states jy^) are collected 
into a matrix T and the measurement vectors \pi) into a 
matrix M, i.e., 

r= [|yo),---,|7m-l)], M =[\po),...,\prn-l)] (7) 

a matrix relation T = MX holds. The entries of the m x 
m matrix X are given by Xki = {p,k\ji)- The joint input- 
output probabilities become p{i,j) = P[S' = i, R = j] = 
and the corresponding correct decision probability 
turns out to be 

m— 1 

Pc=J2 ■ ( 8 ) 

i=0 

The optimization problem can be conveniently formu¬ 
lated introducing the Gram matrix G of the weighted 
states |yi), namely = (yzly^)- Then, using ^ gives 

m—1 m—1 

Gij = ''^^{ji\pk){P'k\jj) = ^ki^kj ( 9 ) 

k=0 k=0 

and G = XtA in matrix form. In conclusion, the op¬ 
timal detection problem reduces to find the factoriza¬ 
tion G = X^ X of the Gram matrix that maximizes 
probability Q. Provided that Xq is the result of the 
optimization, the measurement vectors are obtained by 
M = PAjC^. Note that the independence of the states 


implies that the Gram matrix is positive definite, so that 
X is invertible. Unfortunately the factorization ^ of 
the Gram matrix leads to a set of quadratic equations, 
which can be solved only by numerical programs, at least 
in general. 

A different approach is based on a very popular subop- 
timal measurement, known as square root measurement 
(SRM), which was introduced by Hausladen and Woot- 
ters |S] and thoroughly discussed by Eldar and Forney 
m- In the framework discussed above, the square root 
measurement corresponds to the factorization G = 
of the Gram matrix with X = G^/^ square root of the 
Gram matrix itself. Under our assumptions both the 
matrix G^/^ and its inverse G“^/^ are definite positive. 
The matrix M of the measurement vectors is given by 
M = rG“^/^ and the correct detection probability is the 
sum of the squares of the diagonal entries of G^/^, i.e., 

m— 1 

^[(Gi/%]^ (10) 

i=0 

Even though this measurement is not optimal in gen¬ 
eral, it exhibits interesting properties. It can be straight¬ 
forwardly obtained by the Gram matrix and gives perfor¬ 
mances near to the optimum provided that the weighted 
states |yi) are almost orthogonal. Indeed Hausladen et 
al. na have shown that SRM (also known as ’’pretty 
good” measurement) is asymptotically optimal, in the 
sense that it is good enough to be used as non local mea¬ 
surement in the proof of their fundamental theorem on 
the classical capacity of a quantum channel. More re¬ 
cently SRM has found application HMH in the search 
of optimal measurement for distinguishing hidden sub¬ 
group states, providing an interesting approach to the 
large class of quantum computational problems linked to 
the hidden subgroup problem. In any case SRM furnishes 
a lower bound to the optimal performance and it may be 
assumed as a starting point to find the optimal measure¬ 
ment Xq through unitary transformations. 

The SRM turns out to be the optimal measurement 
when the states exhibit a high degree of symmetry. Then, 
particular attention has been given to applications of the 
SRM to states enjoying geometrically uniform simmetry 
(GUS), i.e., sets of states that are invariant with respect 
to a unitary transformation Hi This is particu¬ 

larly interesting for practical quantum communication 
systems using quantum phase shift- keying (QPSK) |16| 
and pulse position modulation (PPM) [T7]. Also exten¬ 
sions to mixed state have been considered Hi Hi. 

The aim of the present paper is to show that the op¬ 
timal feature of SRM is by no means confined to states 
enjoying GUS and having equal probabilities. In par¬ 
ticular we consider a set of states composed by distinct 
constellations of states enjoying the same GUS (we use 
the term “constellation” borrowed from the telecommu¬ 
nications jargon for the more general “subset”). This is 
a particular case of the compound geometrical uniform 
symmetry discussed in |18| . The novelty of our approach 
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stays in the fact that we find necessary and sufficient con¬ 
ditions in order that SRM is the optimal measurement for 
this case. Moreover we present possible applications to 
practical quantum communication systems as PSK and 
suggest a version of PPM improving its efficiency. 

The paper is organized as follows. In Section II we 
revisit two important results concerning our topic. The 
first, due to Helstrom [50], characterizes the factorization 
G = X that maximixes the correct decision probabil¬ 
ity. The second, due to Sasaki et al. (21], gives simple 
sufficient conditions guaranteeing that the SRM is the op¬ 
timal measurement. In particular the last result enables 
one to show in a very simple way that SRM is optimal 
for states enjoying geometrical uniform symmetry. 
The main result of the paper is presented in Section III, 
where a necessary and sufficient condition is furnished in 
order that the SRM is optimal for a set formed by con¬ 
stellations of GUS states generated by the same unitary 
transformation S but applied to different states. In Sec¬ 
tion IV the result is applied to a double constellation of 
quantum PSK states, showing that SRM may be optimal 
also for states with non uniform probabilities. In Section 
V we consider a double constellation of PPM states and 
show that the SRM is optimal for this case and leads to a 
communications scheme more efficient than the original 
PPM. Some conclusions close the paper. 


II. THE SRM AS OPTIMAL MEASUREMENT 


While Theorem 1 refers to the general case of mixed 
states, a simple characterization of the optimal measure¬ 
ment for independent pure states is given by the following 
theorem. 

Theorem 2. The factorization G = X"^X maximizes the 
correct decision probability if and only if the following 
conditions hold: 

i) for each i and j 

XuX*, = ; ( 11 ) 

ii) the matrix 

y = xxl (12) 

with 


— cll3'§{Woo, . • . ^ X^ — i Yn—i^ (1^) 


is positive definite. 
Hermitian). 


(Note that (11 \ implies that Y is 

□ 


The theorem, not frequently cited in the literature, 
has been proved by Helstrom in a paper |20| concerning 
an iterative search of the optimal measurement. The 
proof of the theorem, given in appendix of the paper, 
is somewhat intricate. A simplified version of the proof 
is given here in Appendix A. 

As we pointed out, the SRM corresponding to the fac¬ 
torization G = X'^X with X = is not optimal in 


general. On the basis of the Theorem 2, Sasaki et al. 
|21| in a paper about the superadditivity of the capacity 
of a quantum channel have found a nice sufficient con¬ 
dition for the optimality of SRM. The following theorem 
gives a generalization of the result. 

Theorem 3. Assume that the Gram matrix G (and its 
square root G^/^) is block diagonal, namely, G = Gi 0 
... 0 Gr. Then, the square root measurement is optimal 
if and only if the square root G^"^ of each block has equal 
diagonal entries. 


Proof. If the matrix G is formed by a single block, 
it is irreducible and its associated directed graph is 
strongly connected [55]. Then, the indexes 0,... ,m — 1 
may be ordered into a cycle iQ,ii, ■ ■ ■, im-i, io such that 
(G^/^)i^i^^^ yf 0 for each k = 0,... ,m — 1. Since G^/^ 
is Hermitian, = (G^/^)and condition 

( flit is satisfied if and only if (G^^^)ikik = (G^''^)ik+iik+i 
for each k or, equivalently, if the diagonal entries of G^G 
have a common value g. Therefore, condition (121 of 
Theorem 2 holds in that XX^ = gG^G is positive def¬ 
inite. If G has many diagonal block, the proof can be 
applied to each block. □ 


Gonfining our attention to the case of a single block, 
since Yu = \Xn\'^ = p(i,i), the theorem is equivalent to 
say that the SRM is optimal if and only if the probabil¬ 
ities of correct decision is g"^, independent of the state 
transmitted. In particular, the optimal correct decision 
probability turns out to be 


Pc = mg"^ . 


(14) 


It is worthwhile to note that the optimality of the SRM 
depends only on the Gram matrix. Since this is invari¬ 
ant with respect to unitary transformations, if a constel¬ 
lation of weighted states satisfy the conditions of The¬ 
orem 3, any other constellation obtained via a unitary 
transformation admits optimal SRM with the same error 
probability. Of course the optimal SRM operators vary 
according to the unitary transformation. 

The most known case of optimality of the SRM con¬ 
cerns states equipped with geometrical uniform symme¬ 
try (GUS) Unilin]. The weighted states ly^) enjoy GUS 
if | 7 i) = S'*| 7 o), i = 0,... ,m — 1, where S' is a unitary 
operator satisfying the condition S'" = I. This general¬ 
izes the usual definition of GUS states jy^) as satisfying 
the condition jy^) = S'|yo), and having equal probabili¬ 
ties Qi = 1/to. The proof of the optimality of the SRM 
is an immediate consequence of Theorem 3. Indeed, the 
entries of the weighted Gram matrix become 

G.j = (7il7i) = (7o|5^-i7o) (15) 

and depend only by (j — i)mod m- One concludes that the 
matrix G is circulant |55j. (Details on the basic proper¬ 
ties of circulant matrices are collected for convenience in 
Appendix B). Then G has spectral decomposition 

G = FAF^ 


(16) 
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) where F is the unitary Fourier matrix (B4) 


(see 

and the diagonal matrix A = diag{Ai,..., Am} collects 
the (positive) eigenvalues of G. Then 




(17) 


is circulant, its diagonal entries are equal, and the SRM 
is optimal. In particular 

1 , m 

g = —Ti (G^/^) = —Tr (A^/^) = — y (18) 

and 

^ ('S') 

The GUS is by no means the only case where Theorem 
3 finds application. As an elementary example consider 
the Gram matrix 




A 

1 


( 20 ) 


of two states | 7 o) and jyi) with equal probability and 
X = (ToIti)- Note that this matrix is not circulant in 
general. Its square root matrix is 


G1/2 


1 

71 


a b 
b* a 


( 21 ) 


with + l&p = 1 and 2ab = x- By virtue of Theorem 3 
the SRM is optimal. A simple algebra shows that 


a=^7l + 7r^, (22) 

so that the Helstrom bound 

Pc = l{l + 7l-|xP) (23) 

holds. This is perhaps the simplest proof of the Helstrom 
bound, at least for equilikely states. In the literature 
the SRM is usually presented in connection with GUS 
states having equal probabilities. However, the exten¬ 
sion to other situations is possible (see for a detailed, 
but scarcely constructive, theoretical analysis). Recently 
Kato pS] has found a closed form expression of the input 
probabilities for which the measurement operators ob¬ 
tained from the not weighted Gram matrix are optimal 
in the case of three coherent states | — a), |0) and |a). 


III. MULTIPLE CONSTELLATIONS OF GUS 
STATES 


Here we consider a collection of sm states {|7fci)}j 
k = = 0,..., TO — 1, defined as 


where S' is a unitary operator such that S'" = I. In 
other words the set of states is formed by s constella¬ 
tions of GUS states obtained by the same unitary op¬ 
erator S, starting from s different states |7fco)- More¬ 
over we relax the hypothesis of equal probabilities, only 
assuming that the states of each constellation have the 
same probability, namely, q^i = Qk with = Ijm- 

Under these assumptions we find a sufficient and neces¬ 
sary condition for the optimality of the SRM. This is a 
particularization of the concept of compound geometrical 
uniform (GGU) states, introduced by Eldar et al. |18j . 
where \xki) = 5'i|7feo) and the to unitary operators Si 
form a (not necessarily Abelian) multiplicative group. 

We begin by ordering the states as 

l7io), ■ • ■ , l7l.m-l), ■ ■ • , |7so), • ■ • , l7s,m-l) • (25) 


The weighted Gram matrix of the states may be parti¬ 
tioned in the form 


G = 


Gii ... Gi, 
Gsi ... Gss 


(26) 


where the submatrices Ghk have order to and entries 


Ghk{i,j) = i^ilXkj) = y/qhQk {ih^lXkj) 

= ^/qhqk{'iho\S^~"\xko) ■ 

It follows that the matrices Ghk are circulant and are 
simultaneously diagonalizable as Ghk = FKhkF\ where 
F is the Fourier matrix defined in Appendix B. Then 
G = FsAFj, with F^ = F 0 ... 0 F and the blocks of 


A = 


All • ■ ■ ^is 
Asi ... Ass 


(27) 


are diagonal matrices of order to. 

The square root of A, say E = A^G^ has the same 
structure as A. Indeed, by a joint permutation of rows 
and columns generated by a permutation matrix H the 
matrix A can be transformed into a block diagonal matrix 


D = HAH^ 


Di ... 0 

0 ... Ds 


(28) 


Incidentally, the permutation converting the ordering 
(251 of the states into the ordering 


l7io),---,l7so),---,l7i ,m—l); • • ■ 1 \^s,m—l} 


(29) 


is the natural choice. Since the matrix A = F^GFs is 
positive definite as G, also D = HAH^ and its diagonal 
blocks Dk are positive definite. Then 


F)G2 = UA^Glf = 



0 

1/2 


(30) 


\lki) = S"\xko) 


(24) 


0 ... D: 
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is well defined. Applying to the inverse permutation 
one gets G^!'^ = where 


Applying the discrete Fourier transform (B5l gives the 
eigenvalues of the block Gn 


b}!"^ = E = 



S 


si 



(31) 


where the blocks T,hk are diagonal matrices of order m. 

The diagonal blocks of the square root of the weighted 
Gram matrix are {G^^‘^)hh = F^hhF^■ This results in a 
circulant matrix with diagonal elements, say gh, depend¬ 
ing on h. The condition of optimality of the SRM (see 
Theorem 3) is satisfied if and only ii gh = g is indepen¬ 
dent of h. This is summarized by the following theorem. 

Theorem 4. The SRM is optimal for the multiple con¬ 
stellation of GUS states if and only if the diagonal blocks 
of the square root G^^^ of the weighted Gram matrix 
have equal diagonal entries. □ 

An alternative criterion derives from the fact that 

gh = -Ti{{G^/^)hh) = -TiiFEhhF^) = -Tr(E;,^) . 
m m m 

(32) 

Then the condition of optimality of the SRM is satisfied 
if and only if the blocks Yihh have equal traces. 

Finally, provided that the SRM is optimal, the correct 
decision probability becomes 

Pc = msg'^ . (33) 


IV. DOUBLE QUANTUM PSK 
CONSTELLATION 

As a first application of the above theory we consider a 
double quantum binary phase shift keying (BPSK) with 
four coherent states | ± a) and | ± /?) in a Fock space. 
The system is a particular case of the scheme discussed 
in the previous section by setting 


Ao = p( 1 + 9a), 

Ai = p[l - ga). 

(38) 

Similarly, the eigenvalues of 

Gi 2 are 


po = Vw(x + 0> 

II 

1 

1 

PX 

(39) 

and the eigenvalues of G 22 



Wo =p(l + 9 / 3 ), 

Wi =p{l-gp). 

(40) 


The weighted Gram matrix is decomposed as G = 
F 2 h.Fl, where A is the matrix ( |27| with diagonal blocks 
All = diag{Ao,Ai}, A 12 = diag{^o,A^i} and A 22 = 
diag{wo, wi}. 

The square root E = A^/^ can be evaluated taking 
advantage of the sparsity of the matrix A, or, as an al¬ 
ternative, considering the permutation Ft that converts 
the ordering of the states into (29). The result is the 


block-diagonal matrix (281, which is related to A through 


D = nAn^, and has diagonal blocks 

Di = 

p(i + 9a) ^/mix + 0 

.a/p9(x + C)* 9 ( 1 + 9 / 3 ). 

D2 = 

p(i-9a) ^/mix-0 

.a/p9(x-C)* 9 ( 1 - 9 / 3 ) _ 


(41) 

(42) 


The square root has the same structure of D, and 
can be evaluated from the square root of and F)\l'^, 
with 


= 

\Jp{^ + Va) + 9(1 + Vp) + 2 V A(+) 

^ p(l -k r]a) + VA(F y/pq{x + 0 _ 

Vmix + O* 9(1 + »?/ 3 ) + Va(+) 

(43) 


I710) = |a), I711) = I - a), I720) = 1 / 3 ), I721) = I - / 3 ) • 

(34) 

In this case the symmetry operator is the rotation oper¬ 
ator S = with a and of annihilation and creation 

operators in the Fock space giving S'!a) = | — a). Accord¬ 
ing to our previous assumption the probabilities of | ± a) 
and I ± P) are p and q respectively, with p + q = 1/2. 
The Gram matrix depends on the inner products 


X = (a|/3), Va = (a| - a), gp = {P\ - P), f = (a| - /?) . 

(35) 

The blocks of the weighted Gram matrix turn out to be 


1 Va 
Va 1 


Gii=p 

Gi 2 = y/pq 


X ? 
i X 


G22 — 9 
= G 


1 9/3 
9/3 1 


21 - 


(36) 

(37) 


where A(+) = pq[[l + 77„)(1 -k gp) - |x -k Similarly 
the square root of D 2 , which can also be obtained from 
(431 with the substitutions 


-9a, 9/3 -t -9/3, C -t -C- 


(44) 


Finally, the square root E = ig obtained applying 
the inverse permutation 11^ to the block-diagonal matrix 
The blocks of E are still diagonal, and in particular 


Ell = diag 


p(l -k ga) + VA(+) 


^p{l + go) + g(l -k gp) + 2 \/aW 

p(1 - ga) + VA(-) 

\Jp{l - 9a) + 9(1 - 9/3) + 2 VA(“) 


( 45 ) 
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and 


E 22 = diag 


g(l + T]/^) + \/A(+) 


p(l + Va) + g(l + Vp) + 2VA(+) 


g(i - ■'7/3) + 


•^^(1 - r]a) + q{l - rjis) + 2 ^^^ 


(46) 


with A^ '> obtained from A*^+) with the substitutions 
(44l. The square root of the weighted Gram matrix is 


then obtained through the Fourier transform, = 

The optimality of SRM can be verified with the aid 
of Theorem 4, or, as an alternative, through the condi¬ 
tions (32). Given a set of parameters a, (3 and p, the 
optimality is not assured. However, it is worthwhile to 
investigate how one of the parameters should be opti¬ 
mized given the others, in order to meet the optimal¬ 
ity conditions. This is the topic of this section, where 
we investigate the optimization of the variable p given 
the parameters a, /3. Although this optimization can be 
(numerically) performed for arbitrary a, j3, we focus our 
attention on two particular cases. 


A. Case |q:| = \j3\ 


As a first case, we consider |Qf| = \I5\. Practical appli¬ 
cation of this case arises when the transmitter, supposed 
to produce a Quaternary-PSK, has a misalignment or a 
systematic bias error in the angle defining one of the two 
constellations. Without loss of generality, the alphabet 
can be defined with two real parameters, a € IR and 
/3 = ae*'^, with 5 € [0,7r/2] defining the angular shift 
between the two binary constellations. 

The entries in the weighted Gram matrix G simplifies, 
and the inner products between the states become 


X = e 




) Va — ■' 7/3 — 


', C = e 




It is immediate to verify that with these positions, con¬ 
dition ( [ 3 ^ for the optimality of the SRM is satisfied by 
p = q = 1/4. That is, the equal prior probability be¬ 
tween the symbols gives the optimality of SRM, even in 
the case 5 f. 

The probability of correct detection is obtained em¬ 
ploying (32) and (33), resulting in 



(\/l + ?7 q + lx + Cl + \/l + 37a — lx + ^1 
+ \/l-3?a + |x-CI + IX-'fl) 


In Fig. [2 the performances of the Double BPSK are 
plotted as a function of the mean photon number |ap 
employed in each transmitted state. In the figure, dif¬ 
ferent curves correspond to the values of the phase shift 
(5 = 0, tt/S, 7r/4, 37r/8,7r/2. As we can see, the lines 



FIG. 1: (Color online) Probability of correct detection for a 
Double BPSK constellation, as a function of the mean photon 
number |a|^ employed in the transmission of each state. From 
the top to the bottom, the performance of different values of 
S as they appear in the legend. 


are monotonically increasing in both the value of \a\’^ 
and 5. In addition, the performances are similar when 
(5 is around 7r/2, which is the situation of the Quater¬ 
nary PSK, while they drop quickly for lower values of 
the phase shift. 


B. Case = 3a 


As a second case of study, we consider o: G M, /3 = 3a. 
This set of constellations corresponds to the modulation 
scheme known as 4-Pulse Amplitude Modulation (PAM). 
While usually an equal prior probability is assumed for 
the transmitted states, here we study the optimization of 
this probability such that the SRM are optimal. 

The inner products between the states read 


X = (a|3a) = 

37/3 = (3a| - 3a) = = rjl, 

C = (a| - 3a) = = pi- 


Again, we resort to conditions (32) in order to find the 
value of p that makes the SRM optimal. In this case, 
there are no evident solutions, and the authors could not 
find any closed form solution, so we turn to numerical 
algorithms to find the optimal value. 

The results of the numerical optimization are plotted 
in Fig. 1^ as a function of the mean photon number |ap 
of the states |±a). Note that while for low values of 
a the prior p moves away from equal distribution, the 
optimal solution tends to p = 0.25 as the value of a 
increases. The corresponding performance is plotted in 
Fig. HI which show an increase in the performance as the 
value of lap increases. 
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FIG. 2: (Color online) Optimal prior probability p of the state 
|±a) of the Double BPSK modulation, in the case a G R, 
P = 3a. 



0 0.4 0.8 1.2 1.6 2 


FIG. 3: (Color online) Probability of error of a 4-PAM system 
with optimized a priori probability and Square Root Measure¬ 
ment at the receiver. On the x axis, the mean photon number 
employed in the transmitted state |±a). 


We have compared the performance of the SRM with 
optimal a priori probabilities with the performances of 
both the SRM and the optimal measurement (evaluated 
via semidefinite programming 0 ) with equal a priori 
probabilities. The results are practically indistinguishible 
if plotted as in Fig. This further comfirms that SRM 
is really a very good measurement in cases of practical 
interest. In particular, in the present case, optimizing 
the a priori probabilities does not reduces in significant 
way the error probability. 


V. DOUBLE QUANTUM PPM 
CONSTELLATION 


A. Simple quantum PPM constellation 


A popular scheme proposed for deep space quantum 
communications is the pulse position modulation (PPM) 
pS] . This scheme can be modeled in a Hilbert space 77®"* 
formed by m replicas of the Fock space H. The states 
used by Alice are m tensor products 

\li) = |7io) 0 ■ . ■ 0 |7i.m-l) (47) 

where only coincides with a coherent state |a), while 
the other coincide with the null state |0). For instance, 
for m = 3, the states are 

l7o) = |a) G) |0) 0 |0), 

|7i) = |0)®|a)®|0), (48) 

I72) = | 0 ) 0 | 0 ) (g) |a). 


In practical realizations, in correspondence with the sym¬ 
bol i the laser radiates only in the 7-th slot of a time 
frame of m slots. 

Although the symmetry of the states is apparent (every 
state is a cyclic permutation of the preceding), the eval¬ 
uation of the unitary symmetry operator S in the tensor 
space 77®"* is by no means trivial [T7]. In any case, pro¬ 
vided that the states have equal probabilities 1/m, the 
weighted Gram matrix is given by 



m 


1 X X 

X 1 X 


X X 1 


(49) 


with X = (a|0)(0|a) = (without loss of generality a 
is assumed real). Since is circulant as G, the SRM 
is optimal. Applying the discrete Fourier transform (B5) 
gives the eigenvalues of G 


Ao = —[1 + (m - l)x], (50) 

m 

Afc = —(I-x), fc = l,...m-l. 

m 


Applying to the eigenvalues of G^/^ the inverse dis¬ 
crete Fourier transform (B6), one gets the first row of 
G1/2, namely. 


Co = 


1/1 -b (m- l)x + (m - 1)^1 - X 
1 


Cl — ... — (^m—1 — /— 

my m l 


- \/r^ 


X 


Note that G^/^ has the same structure of G with equal 
entries out of the diagonal. 

The correct detection probability is obtained by (191 
and is 


Pc = mcQ = 


P _ -i 2 

a/1 -b (m - l)x + (m - l)y/l - X 



















































B. Double quantum PPM constellation 


El = diag{^o,6, • ■ ■ > ^m-i} with 


Here we propose a new PPM scheme with the goal 
of doubling the number of the states and, possibly, of 
improving the capacity of the quantum communication 
system. We consider a double constellation in 

l7lo)j • ■ ■ J l7l,m-l) |72 o)i ■ • • I l72,m-l) ■ 

where, as in the ordinary quantum PPM, \'yii) includes 
only one non null state |a) in the iMh slot, while \^ 2 i) 
includes only one non null state | — a) in the i-th slot. To 
maintain the symmetry of the problem, we assume that 
the states have equal probability 1/(2 to). The weighted 
Gram matrix turns out to to be 

\H K 
K H 


where Hij = (7oi|7oi) = (7ii|7ii) and Kij = (7oi|7ii) = 
( 7 ii| 7 oj). Simple considerations lead to the following cir- 
culant matrices 



'1 

X • 

■ X 


X^ 

X • 

• X 

1 

X 

1 . 

■ X 

, K = — 

X 

X^ • 

• X 

2m 

A 

X • 

. 1 

2m 

.A 

X • 

’• A^ 


The discrete Fourier transform (B51 gives the eigenvalues 
of iJ 


^ = ^[l + (w-l)x], (51) 

fc=l,...,TO-l. 

and the eigenvalues of K 

Mo = + (" 1 -l)x], (52) 

Mfc = ^(x^-X), A: = 1 . 

Then G = F 2 AF'j with 


A = 


Ao Ai 

Ai Aq 


(53) 


where Aq = diagjAo,..., Am-i} and Ai 
diag{/io,..., It follows that 





R T 
T R ’ 


(54) 


where the diagonal matrices Eq and Ei satisfy the con¬ 
ditions Eg -1- Ef = Ag and 2EgEi = Ai. Since the di¬ 
agonal blocks of E coincide, the SRM is optimal. Sim¬ 
ple computations lead to Eg = diagji/Oj i and 


Z/g = 

Vi = 
^0 = 
^1 = 


1 

\/8m - 
1 

\/8m 

1 

VSm 

1 

VSm - 


\/l + X^ -f 2 (to- 1)x+ V^T-^ 


1-X+ V^l-^ 


■\/l + X^ -k2(m- l)x- - X^ 


1-x- 


Finally, using the inverse discrete transform ( |B6[ ), we get 
the first rows of the circulant matrices R and T 


1 


J'o = 


to — 


2my/2m - 


2m\/2m - 


Ti = ti = 


2m\t2m - 


\/l-k2(m-l)x + X^ 

+ {m - 1)(1 - x) + x^ 

■v/l + 2(m- l)x + X^ 

+ im - 1)(1 - x) - x^ 

a/1 + 2(to- 1)x + X^ - (1 - X) 

z = 1,..., m — 1 . 


Note that R and T have equal entries out of the diagonal 
as a consequence of the particular symmetry. 

Finally (331 gives the correct decision probability of 
the double PPM 


Pr = 2mri = 


0 = 2rP [v^l + x" + 2(m-l)x 
-k(m - 1)(1 - x) + m\/l- x^ 


C. Mutual information comparison 


In Fig. I^the error probabilities of the simple and dou¬ 
ble PPM are compared for m = 2 and to = 16 in terms 
of the mean number of photons per symbol related to 
the parameter X by x = s “ . The error probability of 
double PPM is larger, but its states are twice as many. 
Then a significant comparison requires the evaluation of 
the informations transferred by the systems. 

The channel defined by the simple PPM is a symmetric 
channel with joint input-output probabilities p (z, z) = Cg 
and p (z, j) = cf for j ^ i. The marginal probabilities are 
uniform p (i) = p(j) = Ijm. Then the mutual informa¬ 
tion turns out to be 


h 


m—1 

^ p(i,j)\og 
i,j=0 


Pjhj) 
p{t)p U) 


= 2 log TO -I- TOCg log Cg -I- m(m 


(55) 


l)ci log Cl 


The channel defined by the double PPM has joint 
input-output probabilities p{i,i) = Cg, p{i,m + i) = 
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FIG. 4: (Color online) Probability of error of m-PPM 
(dashed) and m-Double PPM (solid), for m = 2 (lower lines) 
and m = 16 (upper lines). On the x axis, the mean photon 
number employed in the transmitted state, |ap. 


FIG. 5: (Color online) Mutual information between the trans¬ 
mitted and estimated symbols, for a m-PPM (dashed) and 
m-Double PPM (solid) constellation, in the case of m = 2 
(lower lines) and m = 16 (upper lines). On the x axis, the 
mean photon number employed in the transmitted state, |a|^. 


p (m + i,i) = Iq, i = 0 ,..., m — 1 , while the other 
4m^ — 4m crossover probabilities have common value . 
The marginal probabilities are uniform p(i) = p{j) = 
l/(2m) as a consequence of the symmetry. In conclusion 
the mutual information turns out to be 


i,j=0 


' p{i)p U) 


l 2 = p{hj) log 2 log(2m) -h 2mrl log 


-I- 2mtQ log tg -I- 4(m — l)mr| log 


(56) 


Asymptotically, as a increases (and y tends to zero), 
Ii tends to logm and I 2 tends to log( 2 m) with a gain 
of one bit per symbol. The mutual informations Ii and 
I 2 , coinciding with the capacities of the channels, are 
compared in Fig. [^for m = 2 and m = 16. 

Two remarks are adequate. First, the double PPM 
scheme is a combination of the simple PPM scheme and 
of a binary PSK scheme. Further gain could be obtained 
combining PPM with fc-PSK, using k PPM constellation 
with non empty states r = 0,1,..., /c— 1. Sec¬ 

ond, from a practical point of view, double PPM implies 
some complications, both on the transmitter and on re¬ 
ceiver side. At the transmitter the on-off modulation 
must be replaced by a phase modulation, while at the 
receiver a simple photon counter must be replaced by a 
phase sensitive device. We do not insist on these topics 
that are beyond the scope of the paper. 


VI. CONCLUSIONS 

The square root measurement furnishes an alterna¬ 
tive approach to the reliable discrimination of quantum 
states. This measurement is only suboptimal, in general. 


but it is well known that it turns to be optimal for quan¬ 
tum states enjoying geometrically uniform simmetry. In 
this paper, we showed that the square root measurement 
is optimal also in other situations of practical interest. In 
particular we found necessary and sufficient conditions 
for the optimality in the presence of multiple constella¬ 
tions of symmetrical states also with non uniform prob¬ 
abilities. As an application example we considered the 
case of two pairs of quantum binary symmetrical states 
(PSK states). Finally, the theory is applied to a possi¬ 
ble improvement of the pulse position modulation (PPM) 
scheme that increases the capacity of the resulting chan¬ 
nel. 
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Appendix A: Proof of Theorem 2 

We apply the Theorem 1 to the operators Wi = | 7 i) ( 7 ^] 
and Pi = \fj,i){p,i\ and preliminarily note the following 
matrix representations of Y, Wr and = Y — Wr in 
terms of the orthonormal basis {|/ij)}: 

{p^\Y\p,) = W,AF , {p^\Wr\pJ) = X,rX*, , 

Y}p = w, AT - X„x*^ . (Al) 

If the measurement is optimal statement i) follows di¬ 
rectly from ([^ and Y = XX^ is semidefinite positive. 
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Since X is not singular, it remains to prove that the di¬ 
agonal entries of Xd (of X) are different from 0. Assume 
that Xqq = (^oIto) = 0- Then, it must exist j ^ 0 such 
that Xj = ^ 0, for otherwise all the states | 7 j) 

would be orthogonal to |^o) against the assumption of 
their linear independence. Without loss of generality, as¬ 
sume that Agi 7 ^ 0. As a consequence of i) we have also 
All = 0. Now, if we change the roles of the measurement 
vectors |/xo) and \fJ-i), the matrix X is modified in A, with 
Xjj = Xjj for j > 1, Ago = Aig, and An = Agi. De¬ 
noted by Pc the new correct decision probability, we get 
Pc — Pc= |Algip -I- |Aigp > 0 contradicting the optimal¬ 
ity of A. One concludes that Xd has non zero diagonal 
entries and AA^ is definite positive. In order to prove 
the sufficiency note that the matrices are Hermi- 
tian and Y = -|- |7r)(7r|- Then they satisfy the 

conditions of a theorem on the eigenvalues of Hermitian 
matrices |221 Theorem 4.3.4] which states that, denoted 
by i^i the eigenvalues of and by Ai the eigenvalues of 
A arranged in increasing order, we have 

1^1 < Ai < < • ■ • < < Am . (A2) 

If A is positive definite, it follows that Ai >0. Since 
in A^’’^ the r-th column vanishes, at least one of its 
eigenvalues vanishes and i/i = 0. One concludes that all 
eigenvalue of A^’'^ are non negative and A^’'^ is semidef- 
inite positive. Since this holds true for each r, the con¬ 
dition of Theorem 1 is satisfied and the operators 
Pi — provide the optimal measurement. 


A matrix G of order m is said circulant if 



Co 

Cl . 

• • ^m— 1 

^m —1 

Cg . 

■ • Cj7i_2 

Cl 

C2 . 

■ • Co 


i.e., if 


G^d = 


2 )mod m 


(Bl) 


(B2) 


In other words G is circulant if its rows are cyclic per¬ 
mutation of the first row. A matrix G is circulant if and 
only if has spectral decomposition 

G = FAA^ (B3) 


where F is the (unitary) Fourier matrix with entries 

= /i,fc = 0,...,m-l (B4) 

\/m 

and the diagonal matrix A = diag{Ag,..., Am-i} collects 
the eigenvalues of G. The circulant matrices of order m 
form a multiplicative commutative group of matrices si¬ 
multaneously diagonalizable. The eigenvalues of the cir¬ 
culant matrix G can be obtained from the first row of the 
matrix via the discrete Fourier transform 


m — 1 

Xk=J2 . (B5) 

r=0 


Appendix B: Circulant matrices 


Finally, the inverse discrete Fourier transform gives the 
first row of G 


We collect in this appendix for convenience some prop¬ 
erties of the circulant matrices which are used in the pa¬ 
per. For more details the reader is deferred to the liter¬ 
ature on the topic, for instance |23| . 


-i m— 1 

c, = - y (gg) 

m 

k=0 
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